1985 


N86 


'24543 

NASA/ASEE SIMMER FACULTY RESEARCH FELLOWSHIP PROGRAM 

Marshall Space Flight Center 
The University of Alabama 


MISSION ANALYSIS FLOW 

Prepared By: 

Academic Rank: 

University and Department: 

NASA/MSFC: 

Laboratory: 

Division: 

Branch: 

Date: 

Contract No: 


SEQUENCING OPTIMIZATION 

Meckinley Scott, Ph.D. 
Professor 

The University of Alabama 
Department of Mathematics 

Program Development, 
Preliminary Design 

Mission Analysis & Integration 

Operations Analysis 

July 12, 1985 

NGT 01-008-021 

The University of Alabama in 
Huntsville 


XXXVI-i 


MISSION ANALYSIS FLOW SEQUENCING OPTIMIZATION 


BY 

Meckinley Scott 
Professor of Mathematics 
Unuversity of Alabama 
Tuscaloosa, AL 

ABSTRACT 

This investigation is an extension of last year's project deal- 
ing with the problem of optimal use of ground resources for future 
space missions. This problem was formulated as a linear programming 
problem using an indirect approach. Instead of minimizing the 
inventory level of needed ground resources, we minimize the over- 
lapping periods during which the same types of resources are used by 
various flights. The model was built upon the assumption that 
during the time interval under consideration, the costs of various 
needed resources remain constant. Under other assumptions concern- 
ing costs of resources, the objective function, in general, assumes 
a non-linear form. In this study, one case where the form of. the 
objective function turns out to be guadratic is considered. Also, 
disadvantages and limitations of the approach used are briefly 
discussed. 
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1. Introduction 

The problem of optimal utilization of ground resources for 
scheduling future space missions has been one of continuing interest 
and concern in the Program Development Division of NASA. Currently, 
planning for efficient use of ground resources is carried out using 
"GROPE" (Ground Resources Operations Program Executive). This is a 
series of computer programs that works with ,a traffic model and 
ground processing time lines as basic inputs. The traffic model is 
the specific Shuttle flights in a given year. Requirements for 
different ground resources and equipment are determined by the 
specific type of each flight, e.g. , a Spacelab (pallet or module), 
deployed satellites with or without upper stages or Department of 
Defense flight. All Shuttle flights require: an Obiter Processing 

Facility where Orbiter refurbishment is done and also where 
horizontally installed payloads are integrated; and a Vertical 
Assembly Building where the solid rockets and external tank are 
stacked on the mobile launcher platform, then the Orbiter is 
attached to this stack, and a launch pad. Many flights also use the 
Vertical Processing Facility for vertically installed payload 
processing. 

GROPE can be constrained or unconstrained . Constrained means 
limited resources are available on limited dates. Here it may not 
be possible for all flights to be scheduled within a given year. If 
this is the case, those flights that are not placed in the schedule 
are moved to the following year and scheduled first. In the 
unconstrained case, the entire traffic model is scheduled in the 
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proper year and it includes a complete specification of the various 
resource requirements in terms of "quantity" and "need" dates. Our 
interest in this project corresponds to the unconstrained case of 
GROPE. A precise description of the actual development of the mech- 
anics of GROPE is, unfortunately, unavailable. However, it is known 
that its development is based mostly on heuristic grounds and lacks 
complete mathematical justification, rigor and formality. Due to 
the presence of a large number of variables, a complete mathematical 
formulation of this problem is no doubt very complex. The problem 
is further complicated by the fact that some of these variables are 
stochastic in nature. In our last year's report [5], we proposed an 
approximate mathematical model to formulate the problem. This form- 
ulation was based on assumptions which describe the actual situation 
fairly closely. It must, however, be pointed out that we have 
ignored a number of variables which have no direct bearing on the 
problem and all variables considered are assumed to be non-stochas- 
tic. Under the assumption that the costs of resources remain 
constant during the period under consideration, using an indirect 
approach, the problem was formulated as a linear programming prob- 
lem. Under other assumptions for the costs of the resources, in 
general, the objective function is non-linear. In this study, one 
case where the form of the objective function turns out to be quad- 
ratic is considered. Also, in general, disadvantages and limita- 
tions of the approach used in this study are discussed. 

2. The Model 

The model was built and dealt with in an indirect manner. 
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Instead of analyzing the problem by minimizing the inventory level 
of needed resources, we use an approach whereby the flights are 
scheduled in such a way so as to minimize the overlapping periods 
during which the same types of resources are used by the various 
flights. Associated with each overlap for the use of the same type 
of resource by any pair of flights is a penalty cost which depends 
upon the number of units of that resource needed by these flights. 
Assuming that costs of resources remain unchanged during the whole 
period of the traffic model under consideration and that penalty 
costs are directly proportional to the lengths of the corresponding 
overlaps, the over-all objective function, which is the sum of such 
costs, is linear. Associated with the objective function is a 
number of sequencing and resource constraints which are also linear 
in form. Thus the problem under consideration falls within the 
domain of linear programming. 

One program involvig use of linear programming for scheduling 
flights concerning space mission was prepared by Lockheed Electron- 
ics Co. in 1976 in the form of a technical report [1], However, use 
of this program is basically restricted to cases dealing only with a 
fixed set of resources and the problem is one of selecting a traffic 
model from among various flight candidates which satisfy certain 
objectives. In our investigation, we have elaborated and refined 
the approach proposed in [6]. Stated below are the basic assump- 
tions that concern our investigation. 

(1) There are n flights to be launched during some given interval 
of time, [0, T]. 
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(2) The order in which these flights are to be launched is prede- 
termined. 

(3) Of the n flights, there are p (-^n) specific flights with fixed 
launch dates. Each such flight may, however, have a launch window 
of a certain specified length. 

(4) All flights utilize at least one type of ground resource from a 
collection of M different types. 

(5) Any flight that utilizes say, a type 'k' ground resource may 
require n^ (^l) units of that resource. 

3. Notation and Formulation 

Let, t^, (i = l,2,...,n) be a variable denoting the launch 
time of the ith flight. 

For any flight i that utilizes a type 'k' resource before its 
launch time, sM denotes the length of time in which this 
resource is seized before t^ and d£l denotes the corresponding 
duration for its use. 

For 1 ^ i < j < n, the non-negative overlap^ variables are 

denoted by 0.K1, (k = 1,2, ...,M; r = 1,2 3,4). 

The variable 0?^ measures the amount of the overlapping 

period for use of a type 'k' resource by the ith and jth flights 

when the seize times for the resource occur before their launch 

times. The variable 0^-2 measures the amount of the overlapping 

period for use of a type 'k' resource by the ith and jth flights 

when the seize times for the resource occur before tf for the ith 

flight and after t-i for the jth flight. The variables O.M and 

J L J 

0^4 are similarly interpreted. In the case of oM, the 
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seize times for the resource occurs after ti for the ith flight 
and before tj for the jth flight. For oM, the seize times 
for the resource by both flights occur after their launch times. 

Clearly, for a traffic model consisting of n flights and M types 
of resources, the launch time variables t^,t 2 f...,t n generate 
2Mn(n-l) non-negative overlap variables. 

The assumption that no two flights are to be scheduled at the 
same instant of time leads to constraints of the type 

ti+l - tj_ ^ d > 0, (i = 1,2, ,n ) . (1) 

The provision that there are p specific flights with fixed 

launch dates at times, say, t n . , (i = 1,2,. ..,n) may have launch 
windows yields constraints of the type 

t n - 1. i f n 1 i # (i — l,2,...,p). (2) 

Also note that t n ^ T. (3) 

The relationships between the variables oH r and t^, 

(1 i <■ j 4 n, r = 1,2, 3, 4; k = 1,2,.. . ,M ) are given by the 
equality constraints 

- ti + tj + off = dK r , (4) 

where the constants dK r are defined by 
dK.l = d^l - s^l + s-kl 

C t J ' 

dK2 = d^l — ski - s^2 

tj t «- u ' 

dK3 = d^2 + s^2 + s^l 

tj t *. j ' 

<jk4 _ <jk2 + s k2 _ s k2 > 

Cj i i J 

Let nkl and n^2 denote, respectively, the number of 

l <► 

units of a type 'k' resource utilized by the ith flight when the 
seize times for this resource occur before and after t£, n^l 
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and n^ 2 for the jth flight are similarly defined. These quanti- 

J 

ties are used in the formulation of the cost functions. Under the 
assumption that costs of resources remain constant during the period 
under consideration, we consider the case where the cost, c£. r 
associated with the overlap variable oH r is given by 


II 

m y 
u 

n kr 0 kr 

tj ij 

C< k >, 

(1 ^ 

i < j n 

; r = 

Here, 

n kl = 

Min(n k l, 

n k l) 

L J 

L 

J 

n k ,2 = 

Min(n k l, 

n k 2 ) 


cj L J 

and C*k) is the weight given to a type 'k 1 resource based on 
cost considerations relative to other types of resources in the 
collection. The over-all objective function T, say, is the sum of 
costs C kr over all values of k, r, i and j, (i 4 j). 

After a considerable amount of straight forward algebra, it can be 
shown that mininizing T is equivalent to minimizing T^ given by 




Ti = 


where 


Ci = 


m 



Citi, 


L - \ 


M 


1 

4 


(k) 







( i 1/ 2^ ( • *.fii ) 


£ = n kr , (1 ^ i < j i 

^ U 

T - » 

Recall that T^ which is linear in t^'s is to be minimized 
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subject to linear constraints given by (1), (2), (3), and (4). 

4. Some Remarks 

In the model developed above, we have assumed that there are p 
specific flights out of n with launch dates fixed at times t n , , 

L 

(i = 1, 2, 3, . . . ,p ) ; each of which has a launch window of a certain 
specified length. Since the order of the flights for any given 
traffic model is predetermined, the launch times t£, 

(i = 1,2,... ,n) satisfy 

0 ^ tx t2 «... <t n ^ ^ ^ 

fc n x < tn a -H » < fc nj + i 

... < t n .j <t n . 

In the case where np = n, that is, the last flight is one with a 
fixed launch date, the constraint t n ^ T will be replaced by one 
of the form 
fc n ^ Ip* 

Since all flights are to take place during the interval [0,T], I p 
should satisfy the condition 

j-p — x . 

Similarly, in the case where n^ = 1, we have the condition 

fc l il 
with li 0. 

5. Other Cost Functions 

In the model developed above, due to the assumption that costs 

of various resources remain constant, the objective function turned 

out to be linear. If this assumption is relaxed, in general, the 

objective function is non-linear and this would considerably 

increase the degree of difficulty of the problem. Below we consider 

the case where the objective function assumes a quadratic form. 
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Here, we assume that the cost of each type of resource increases 
linearly with increasing time. Specifically, we assume that the 
cost C (k) for the kth resource at time t£ is 


C (k) = a k + b k ti. 


and the penalty cost C^f associated with the overlap variable 

0 kr is given by 
l j 

C kr = (C< k > - C< k) ) n kr oH.r 
j «■ ij t -j 

= b k (tj - ti) n kr 0 kr , 

Lj L j 


(r = 1,2, 3,4; k = 1,2,.. .,M; 1 i <- j t^n). 
The overall objective function is thus 


T = V 1 C kr , 

— i L J 

where the sum is to be performed over all values of r, k, i and j, 
(1 ^ i a j < n ) . 


After a lengthy but straight forward algebra, it can be shown that 


T = 



!ij (t i - tj> 2 + X/ij 


(ti - tj) 


where each sum is to be performed over the set 
^(i,j)/l i i i j ^ n J , and 

M q 


13 


£ X s 


*>--=■ « t=. i 

M H 


ij = - V V n k ,r d k ? b k . 

x— — < U C J 


fc-i t--i 
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Further simplification reduces the objective function to the form 


T = ( Pjtj + 2(ljjtjtj) + j^liit 2 , 

bl I £ t < j £ V\ ^ = ' 

where 



• -tci i <-r\ 


The form of the objective function T does not in any way alter 
the form of the constraints and hence T is to be minimized subject 
to the same constraints given by (1), (2), (3), and (4). In terms 
of matrix notation, the objective function can be written as 
T = pt N + tLt v , 
where 

L = (lij) is a symmetric n X n matrix, 
h ( t j,t 2* . . . ft n ) , 

£ = ( PI rP2 • • • • fPn ) • 

and t N and p' are column vectors corresponding to t and £. 

6. Discussion and Conclusion 

At the present time, there are some convincing arguments that if 
we are dealing with a relatively short interval of time, (e.g., one 
year), the costs of resources could very well remain constant. 

Thus, in such a case, the formulation of the problem as a linear 
programming problem is well justified. Furthermore, this assumption 
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does lend itself to easier treatment both from the mathematical and 
computational aspects of the problem. Other assumptions for the 
cost function, in general, considerably increase the complexity of 
the problem especially from the computational point of view. 

As mentioned earlier, for a traffic model with n flights, in 
addition to the n launch times, variables, in general, there are 
2Mn(n-l) overlap variables and the number of constraints is n + 2p + 
2Mn(n-l). Thus, if there are no flights with fixed launch dates, 
the number of variables is the same as the number of constraints. 

An increase in the number of resources by one may cause an increase 
of as many as 2n(n-l) variables. In actual practice, the number of 
variables actually needed may be reduced if, for example, we know 
that the duration for use of any resource is short relative to the 
interval under consideration. On the other hand, for any flight 
type, there would be a minimum of about 15 ground resources where 
cost considerations should not be ignored. Thus, even if we have as 
few as 20 flights per year, the number of variables involved may be 
as many as 20 + 2 (15 ) ( 20 ) (19 ) = 11420, a very large number indeed. 
The recent findings of Karmarkar [2] concerning an algorithm for 
large scale linear programming problems when available will, hope- 
fully, prove useful to our problem. So far his approach exists only 
in what has been described as rougher computer code and, perhaps, it 
has not been tested on a : wide range of problems. 

In the analysis of the problem, we have used an indirect 
approach to optimize the use of needed ground resources. This 
approach provides us with launch times for scheduling flights of any 
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given traffic model. However, note that it does not directly lend 
itself to the computation of the total cost o'f resources needed, 
order to find the total cost corresponding : t'o any 'solution 
( ti,t 2t • • • # t n > °f the launch times, the peak requirement for 
each resource will have to be determined -in each case. 
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